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Abstract: With an effective chiral flavour SU(3) model we show the effect of hadronic resonances on the QCD 

phase diagram. We state that varying the resonance couplings to the scalar and vector fields affects the 
order and location of the phase transition, the possible existence of a critical end point (CEP), and the 
thermodynamic properties. We present (strange) quark number susceptibilities at zero baryochemical 
potential and at three different points at the phase transition. Comparing results to lattice QCD, we state 
that reasonable large vector couplings limit the phase transition to a smooth crossover ruling out a CEP. 
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1. Introduction 

Both on the theoretical and on the experimental side there are great efforts to study the largely unknown properties 
of the QCD phase diagram with special attention to the phase transitions of strongly interacting matter. Most 
interesting are the transitions that restore chiral symmetry and that generate deconfinement from hadrons to a 
Quark Gluon Plasma (QGP) at high temperatures and densities. Experimentally different regions of the phase 
diagram are studied in ultra-relativistic heavy ion collisions at different beam energies. Since QCD cannot be 
solved analytically, there are mainly two groups of theoretical approaches to QCD. On the one hand there are 
lattice QCD models, which are mostly limited to vanishing baryochemical potentials, and on the other there are 
effective models describing specific properties of QCD matter. Here, we present an effective chiral flavour SU(3) 
model with hadronic degrees of freedom including heavy resonance states and study the properties of the phase 
diagram compared to lattice QCD results. 
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Figure 1. (Color online) Left: Normalised chiral condensate order parameter together with lattice data for different actions. 

Right: Phase transition lines for different r v . Dashed lines depict a crossover, solid lines a first order phase transition. 



2. Model 

Our S't7(3)-flavour a-oj-model with a non-linear realization of chiral symmetry [1] is based on the mean field 
Lagrangian C = £kin + £int + £mes- With the hadrons' kinetic energy in the first term, the interaction of baryons 
with scalar a, £ and vector mesons u), (f> is described by 

Cim = - ^2 ipi (m* + gtuijooj + 3i07oc/>°) V»- (1) 

i 

The index i includes the baryon octet, decuplet, and heavier resonances with masses m < 2.6 GeV with a 
minimum three star rating in the Particle Data Book. The term £ mcs = £ V cc + £n + £esb includes the vector 
and the scalar mesons' self interactions and the explicit chiral symmetry breaking. The effective baryon masses 
m i — 9ia a + ftcC + 3 m i are generated by the coupling of baryons to the scalar meson fields a, ( and by an 
explicit mass Sim. At high temperatures and baryonic densities this formalism leads to smaller baryon masses 
and thereby to the restoration of chiral symmetry due to the decreasing <r-field. The couplings of the baryon octet 
to the mesonic fields and the mesonic potential are fixed to reproduce the tabulated vacuum masses and well- 
known nuclear ground state properties. The coupling strengths of the baryonic resonances are set proportional 
to the nucleon couplings gjv via the parameters r s , v according to gBa,u = r s,v • gNo.u and (/sc^ = r s,v • gNC,4>- 
For simplicity and for generating the effective baryon masses dynamically by the scalar fields, we keep the scalar 
coupling fixed at r 3 = 0.97. Thereby, we are able to reproduce the particles' vacuum masses, except for the 
explicit mass 5rrn, and to ensure a smooth crossover transition at zero baryochemical potential [ib = 0. However, 
the vector coupling parameter r v , which controls the abundances of the baryonic resonances at /is , is varied 
in order to study its impact on the calculated phase diagram and the thermodynamic properties. In contrast 
to this interacting hadron resonance gas (HRG), we also perform calculations for an ideal non- interacting HRG. 
In this particular case r 3 ,„ are set to zero and the masses of all particles are fixed at their tabulated vacuum 
expectation value. 
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Fi gure 2. (Color online) Interaction measure (e — 3p)/T 4 (left) and the strange susceptibility Xs/T 2 (right) for the interacting 
(blue lines) and non-interacting (black lines) HRG compared to lattice data. 



The grand canonical potential Q/V = — Ant— £ m eso;n+fith includes the thermal contribution of all hadrons (i.e. all 
particles and antiparticles) in the model. The effective baryochemical potential is defined as [i* = fj,i — gi U ui — gi^cj). 
From the grand canonical potential all thermodynamic quantities follow, i.e. the pressure p, the energy and entropy 
densities e, s, and the particle densities pi. 

3. Results 

The normalised chiral order parameter a/oo(T) from our model is shown in Fig. 1 (left) together with lattice 
data from different actions and lattice spacings [2]. With r s fixed, a/ao of the interacting HRG (blue dotted line) 
exhibits a smooth crossover and a critical temperature T c = 164 MeV being in line with recent lattice data for 
the continuum extrapolated HISQ and stout actions. Analysing a/ao in the whole T — /x-plane we get the phase 
transition lines shown in Fig. 1 (right) for different values of r v . Here, dashed lines depict crossover transitions, 
dots a CEP, and full lines first order phase transitions. We find that for reasonable values of r v > 0.4, for which 
the nuclear ground state is correctly located in the chirally broken phase, the CEP moves to low temperatures or 
vanishes completely (r v > 0.6) in favour of a broad crossover transition. 

Figure 2 shows the interaction measure (e — 3p)/T 4 (left) and the strange susceptibility Xe/T 2 (right) reflecting 
strange quark fluctuations [3] both as functions of T. As expected, compared to the non-interacting HRG the 
interacting HRG (blue lines) shows a much faster increase of both quantities at T c due to more degrees of freedom. 
Both quantities from the interacting HRG are in qualitatively good agreement with the stout continuum limit. 
Studying the second and fourth order quark number susceptibility coefficients C2, C4 at \ib — (Fig. 3 (a) and 
(b)), we find that fluctuations at T c get massively suppressed by large couplings to the repulsive vector field 
r v . This is also reflected in the susceptibility ratio C4/C2 ((Fig. 3 (c)) which only comes close to lattice data for 
reasonably high values of r v > 0.8. This again supports our previously mentioned preference for larger values r v 
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Figure 3. (Color online) Susceptibility coefficients C2, C4 ((a) and (b)) and the ratio C4/C2 as a function of T/T c at fis — 
(c) for the interacting (blue lines) and non-interacting (black lines) HRG for r v — 0.8 and at different points (CEP - 
liB = 216 MeV; CO - £1 B = 30 MeV; PT - /i B = 490 MeV) on the phase transition for r„ = 0.9 with lattice data (d). 



and thus limits the phase transition to a broad crossover in the whole T — ^t-range (cf. Fig. 1 (right)). For the quark 
number susceptibilities we think, a comparison to not yet available lattice data with continuum extrapolated stout 
action would be very promising. A strong suppression of the susceptibility coefficients is also seen at /jb / 0. 
When calculating the ratio C4/C2 crossing the phase transition at different values of [is (Fig. 3 (d)), we find that 
fluctuations are much smaller compared to /is = which is again due to the repulsive vector interactions. 
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